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Abstract: A moving pellicle superconducting mirror, which is driven by radiation pressure 
on its one side, and by the Coulomb force on its other side, can become a parametric oscillator 
that can generate microwaves when placed within a high-Q superconducting Fabry-Perot 
resonator system. A paraxial- wave analysis shows that the fundamental resonator eigenmode 
needed for parametric oscillation is the TM ii mode. A double Fabry-Perot structure is 
introduced to resonate the pump and the idler modes, but to reject the parasitic anti-Stokes 
mode. The threshold for oscillation is estimated based on the radiation-pressure coupling 
of the pump to the signal and idler modes, and indicates that the experiment is feasible to 
perform. 

I. INTRODUCTION 

Parametric oscillators for generating electromagnetic microwaves might be possible, based upon 
the idea that a moving mirror is like a moving piston that can perform work nonadiabatically on 
radiation contained within a cavity. Above a certain threshold, the parametric action of the moving 
mirror will exponentially amplify this radiation until it can become a large-amplitude, classical 
wave. An essential component for such a parametric oscillator is a high-Q superconducting (SC) 
resonator, which can be formed by a fixed, spherical curved SC mirror in conjunction with a moving, 
fiat SC mirror (see Figure 1). 

To understand how the proposed parametric oscillator works, let us imagine performing a 
thought experiment in which some weak, on-resonance, "seed" microwave radiation from an outside 
microwave source, is injected backwards through the output hole of the curved SC mirror of Figure 
1, into the volume of the cavity between the curved and fiat mirrors. Thus the SC resonator can 
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FIG. 1: Hemiconfocal Fabry-Perot superconducting (SC) resonator with electrostatically charged mirrors, 
for generating microwaves in a parametric oscillator. A Gaussian beam mode focused by a curved mirror onto 
a flat mirror with charge q that can move back and forth parallel to the z axis, yields a longitudinal electric 
field E z satisfying conducting boundary conditions. The lowest mode of the resonator is the transverse 
magnetic TMon mode (see Appendix A). 

be filled with "seed" radiation, exciting it into its fundamental TMon eigenmode. As a result, 
there will appear a longitudinal electric field E 2 at the surface of the SC flat mirror that will be 
oscillating at the TMon eigenmode frequency, i.e., at a microwave frequency (see Appendix A). 

Next, imagine that the flat SC mirror consists of a thin SC film sputtered onto the left side of 
a thin, light, flexible diaphragm (e.g., the pellicle mirror to be introduced later in Figure 4), which 
is sufficiently thin so that this diaphragm can easily be driven into mechanical motion. Moreover, 
imagine that this film is electrostatically charged with a DC charge q. Then the longitudinal electric 
field E z at the surface of the SC film will lead to a force 

F z (t) = qE z (t) (1) 

oscillating at the same microwave frequency as that of E z (t) . Note that ([I]) is a linear relationship 
between the force (t) and the electric field E 2 (t). In this way, a force oscillating at a microwave 
frequency can be exerted upon the diaphragm parallel to the z axis, and the flat SC mirror will 
be driven into simple harmonic motion at the same microwave frequency as that of the TMon 
eigenmode of the SC resonator of Figure 1. 

Now consider the opto-mechanical configuration sketched in Figure 2, in which a laser beam 
from the right is incident on a moving mirror (e.g., on the flat SC mirror of Figure 1 with mul- 
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tilayer dielectric optical coatings deposited on its right side). When this moving optical mirror 
is combined with a fixed optical mirror in an optical Fabry-Perot-cavity configuration, there will 
arise a production of Doppler sidebands, which can then be utilized either for the laser cooling of 
the moving mirror by means of a red-detuned laser tuned to the lower Doppler sideband PQ, or for 
the parametric oscillation of the moving mirror excited in an elastic mode at acoustical frequencies 
by means of a blue-detuned laser tuned to the upper Doppler sideband [2]. Above the threshold 
for parametric oscillation of the moving mirror within the SC resonator of Figure 1, a "signal" 
wave would begin to build up, growing exponentially with time starting from the injected "seed" 
microwave radiation. However, once parametric oscillation above threshold has occurred, one could 
turn off the source of the "seed" radiation. The SC resonator would then continue to oscillate as an 
autonomous source of the same microwaves as the "seed" in the absence of the "seed", just like the 
autonomous generation of microwave radiation by the original ammonia maser above its oscillation 
threshold [3]. 
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FIG. 2: A mass-and-spring model for an opto-mechanical parametric oscillator that could generate microwave 
radiation starting from "seed" radiation. The SC Fabry-Perot microwave resonator in Figure 1 is modeled 
by a charged mass attached via a spring to a fixed, charged wall. The "pump" laser beam (in blue) is blue- 
detuned so as to coincide with the upper Doppler sideband arising from the simple harmonic motion of the 
moving mirror. After reflection from the moving mirror, the laser beam is red-shifted to become the "idler" 
beam (in red), whose frequency coincides with the lower Doppler sideband. The "pump" laser beam will 
exert a radiation pressure force that will drive the mirror into parametric oscillations. Here the "Double" 
optical Fabry-Perot resonator represents a simplified model for the doubly-resonant structure depicted in 
Figure 3, which serves to resonate both the "idler" and the "pump" frequencies, but to reject the any Doppler 
up-shifted frequencies. 
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II. OPTO-MECHANICAL MODEL FOR PARAMETRIC OSCILLATORS 

For ease of understanding, consider the simplified opto-mechanical model for the parametric 
process shown in Figure 2, in which a moving mirror attached to a spring is coupled via radiation 
pressure to a strong "pump" laser at frequency u p , a weak "idler" light wave within an optical 
resonator at frequency Wj, and a weak "signal" wave in the SC resonator at frequency oj s , which is 
being represented as a simple harmonic oscillator with a resonance frequency uj s , such that 

u}p = uji+uj s (2) 

The energy for the oscillating fields at u>i and ui s comes from the energy supplied by the pump laser 
at ojp. 

In Figure 2, we have replaced the hemiconfocal microwave resonator shown in Figure 1 by a 
mass-and-spring model, in which the simple harmonic mechanical motion of the flat SC mirror of 
Figure 1 driven in the presence of the charge q by the longitudinal electric field E 2 (t) parallel to 
the z axis, is modeled by the simple harmonic motion of a mass attached via a spring to the fixed 
wall on the left side of Figure 2. The justification for using this mechanical model (see Appendix B) 
is that we have found that the fundamental eigenmode solution of the SC hemiconfocal microwave 
resonator, which is the Gaussian-beam TMon mode, possesses a longitudinal component E z of 
the electric field which has a nonvanishing component along the z axis of the resonator acting on 
the electrostatic charge q at the surface of the flat mirror of Figure 1. Therefore there exists a 
longitudinal component of the force that drives the flat mirror back and forth into simple harmonic 
motion along the z axis, whose motion can then be simulated by the mass-and-spring model of 
Figure 2. 

Furthermore, we shall model, by means of a single, simplified Fabry-Perot structure on the 
right side of Figure 2, a "Double" optical Fabry-Perot resonator, which is depicted in more detail 
in Figure 3, whose purpose is to simultaneously resonate both the strong, incoming blue-detuned 
laser (i.e., the "pump" laser) and the weak, Doppler red-shifted "idler" light wave produced upon 
reflection from the moving mirror, but whose purpose is also to serve as a rejection filter to reject 
any undesirable Doppler blue-shifted (or "anti-Stokes") light. Later, we shall put all these pieces 
of the parametric oscillator system together in Figure 4. 
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FIG. 3: (a) "Double" Fabry-Perot is a three-mirror resonator system with three equally spaced mirrors Ml, 
M2, and M3. The leftmost mirror Ml is the flat, moving mirror attached to the spring of Figure 2, i.e., 
the "pellicle-mirror" of Figure 4. The pump laser is incident from the right through the rightmost mirror 
M3. Multi-layer dielectric coatings are indicated by red lines; anti-reflection coatings by blue lines. The 
reflectivity of M2 is less than that of Ml and M3. (b) The spectral transmission of this resonator system 
has a doublet structure, whose splitting is proportional to the transmittivity of the middle mirror M2. The 
frequency of the "pump" laser of the parametric oscillator coincides with the upper member of the doublet, 
and the frequency of the "idler" ("Stokes") wave coincides with the lower member. The doublet splitting 
frequency is resonant with the "signal" frequency of the parametric oscillator. The parasitic "Anti-Stokes" 
wave is rejected by the structure. 



III. THE "DOUBLE" FABRY-PEROT RESONATOR 



Figure 3(a) shows the detailed structure of the "Double" Fabry-Perot resonator, which consists 
of three equally spaced mirrors Ml, M2, and M3. Suppose that Ml is the moving mirror of Figure 
2 (or the "pellicle mirror" of Figure 4) . Suppose further that at t = light initially fills the left half 
of the structure between Ml and M2, but that there is initially no light in the right half between 
M2 and M3. Then after a period of time determined by the transmittivity of the middle mirror M2, 
light will leak from the left cavity to the right cavity, until light fills the right cavity. This results 
in a periodic sloshing back and forth of the light between the left and right halves of the structure 
at a frequency equal to the splitting frequency of the spectral doublet shown in Figure 3(b). The 
"pump" laser will be tuned into the resonance with the upper member of the spectral doublet; the 
"idler" wave to the lower member of this doublet. Any small amount of noise at the idler frequency 
will lead to a radiation pressure force on the moving mirror Ml that will be modulated at the 
doublet splitting frequency, which can be chosen to be tuned into resonance with the SC microwave 
resonator frequency of Figure 1. This then can lead to a mutual, resonant reinforcement of the 
noise in the "idler" and "signal" modes that leads to exponential growth of both by parametric 
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amplification, and possibly to oscillation above a certain threshold, as will be presently shown. 

An important feature of the "Double" Fabry-Perot is that it automatically rejects all Doppler 
up-shifted, or "anti-Stokes," spectral components arising from the motion of the moving mirror 
towards the laser in Figure 2. Such undesirable "anti-Stokes," or up-shifted, frequency components 
of the light can rob energy away from the desired "Stokes," down-shifted, frequency components, 
specifically, the "idler" frequency necessary for exponential parametric gain to occur. Since the 
spectral doublet structure of Figure 3(b) does not have any resonance at the "anti-Stokes" frequency, 
the "Double" Fabry-Perot will serve as a spectral rejection filter that prevents light from building 
up inside the resonator at any unwanted, Doppler up-shifted frequencies. 



IV. PELLICLE-MIRROR FABRY-PEROT SYSTEM 
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FIG. 4: Pellicle-mirror parametric oscillator for generating microwaves. A thin, charged SC niobium film 
(in black) is sputtered onto the left side of a thin pellicle substrate (in green). A thin, oppositely charged 
SC niobium film (also in black) is sputtered onto the reflecting surface of the leftmost curved SC mirror, in 
order to form an extremely high-Q SC resonator. Optical coatings (in red) are evaporated onto the right side 
of the pellicle mirror. Thus the pellicle mirror divides the entire structure into two halves. In the left half, 
the microwave-frequency Fabry-Perot resonator is for resonating a microwave "signal" within a Gaussian 
beam (in dotted black). In the right half, the "Double" optical Fabry-Perot resonator is for resonating both 
the "idler" light beam (in dashed red) and the blue-detuned "pump" laser beam (in dashed blue). The 
microwave mode on the left side of the pellicle is the TMon mode, and the optical mode on the right side 
of the pellicle is the TEM 00n mode. 



A practical implementation of the "moving mirror" idea at the heart of the microwave parametric 
oscillator, is to evaporate a thin SC film onto a thin, flexible, low-mass pellicle (with a thickness of 
around two microns [4]) stretched tautly over a circular wire frame in order to form a "drumhead" 
with a low-frequency acoustical eigenmode whose mode pattern has a central maximum at the 
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center of the pellicle. On the reverse side of the pellicle, one could then evaporate dielectric optical 
coatings to reflect both the strong "pump" laser light beam, and the weak "idler" light beam, for 
implementing the "Double" Fabry-Perot cavity for the "pump" and "idler" optical beams, which 
is sketched in Figure 3. 

In Figure 4 (which is not drawn to scale, since the diffraction of the microwave- frequency "signal" 
Gaussian beam will occur much more quickly than the diffraction of the optical-frequency pump 
and idler light beams), we put together all the pieces of the parametric oscillator system based on 
the "pellicle mirror" idea. Since microwave frequencies are very much higher than the acoustical 
frequencies of the pellicle drumhead, the motion of the central portion of the pellicle is essentially 
like that of a free body (i.e., a free mass) being driven at microwave frequencies. 

One can understand the parametric amplification arising in this configuration as follows: When 
the "pump" laser beam enters from the right into the "Double" optical Fabry-Perot system, the 
pellicle will act as if it were a piston that is being driven by the radiation pressure force which 
is varying at the beat note frequency between the strong "pump" laser beam (in blue), whose 
frequency is u p , and the weak "idler" (or "Stokes") beam (in red), whose frequency is Wj, such that 
this beat note frequency equals the "signal" frequency u s , which is the TMon eigenmode resonance 
frequency of the left-hand-side microwave Fabry-Perot resonator. Then during the parametric 
amplification process, one quantum of the pump wave will break up into one quantum of the signal 
wave plus one quantum of the idler of wave, satisfying the energy conservation relationship 

hujp = huji + huj s (3) 

where ftw p is one quantum of energy of the "pump" mode, hu>i is one quantum of energy of the 
"idler" mode, and hu s is one quantum of energy of the "signal" mode. 

The moving SC pellicle mirror, when it is viewed as if it were a moving piston, will do work on 
the "seed" microwaves contained inside the left Fabry-Perot as it moves nonadiabatically [5]. Hence 
the action of the moving pellicle will amplify this "seed" radiation, which will in turn amplify the 
motion of the pellicle when it is viewed as a moving optical mirror, because this motion will further 
amplify the strength of the "Stokes" Doppler sideband of the incoming "pump" laser beam, thus 
amplifying the strength of the "idler" beam. This in turn increases the strength of the beat note 
in the radiation pressure force acting on the pellicle, which further increases the amplitude of the 
motion of the pellicle, etc., in a feedback process. In this way, there will be a mutual reinforcement 
of the "signal" wave and the "idler" wave, so that an exponential growth of both waves (above a 
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certain threshold of oscillation) will result, as in a traditional laser. For each quantum produced in 
the "idler" mode, a quantum of the "signal" mode will be produced as well, in agreement with the 
Manley-Rowe relations [6]. 

V. THRESHOLD FOR PARAMETRIC AMPLIFICATION 

We make a simple argument to estimate the threshold for parametric amplification of the pellicle 
mirror parametric oscillator shown in Figure 4. Let us assume that the microwave cavity on the left 
side of the parametric amplifier can be modeled as a mirror with a spring attached to a fixed wall 
(see Figure 2) so that it forms a simple harmonic oscillator with resonant frequency cj s , effective 
mass m, and quality factor Q s . The right half of the parametric oscillator is regarded as a "Double" 
optical Fabry-Perot resonator with two resonances at and uj p . The transmission of this "Double" 
Fabry-Perot has been studied in [8] and is illustrated in Figure 5. Note that the splitting between 
the double peaks depends on the values of the reflection coefficients. We require that n = = r 
and T2 < r, where r corresponds to reflection coefficient of the two end mirrors, and r<i to the 
middle mirror in Figure 3(a). Furthermore, the purpose of the "Double" Fabry-Perot is to allow 
for the selection of the two desired optical modes uji and co p , as already discussed. Since the 
"Double" Fabry-Perot acts as a single Fabry-Perot with two closely spaced resonances, we treat 
the parametric amplifier as a single Fabry-Perot cavity with a harmonically moving end mirror as 
illustrated in Figure 2. 
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FIG. 5: The transmission of a "Double" Fabry-Perot has resonances at the "idler" mode w$ and the "pump" 
mode ujp. The difference frequency between "pump" and "idler" is resonant with the "signal" frequency 
(i.e., io p — u)\ = uj s ). The "anti-Stokes" cj£ is off resonance and, hence suppressed. 

The Fabry-Perot with a harmonically moving end-mirror (see Appendix C) yields the following 
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equations of motion 

^ + 7i £bi = ^X*f 0pe -^ (5) 

where Aw = lo p — uii — u B , A is the cross sectional area of the resonator, X is a slowly varying 
oscillator amplitude, r is the round trip time between the two mirrors of Figure 2, <?oi is the slowly 
varying amplitude of the "idler" mode, and 7i , 7 S are the HWHM for the "idler" and "signal", 
respectively. The problem of a Fabry-Perot with a harmonically moving end-mirror has been solved 
in detail by [9]. Choosing the same normalization (see Appendix C) as in |9], we arrive at the same 
results 

§+7 s X = ^D p D*e-^ (6) 
at muj s L 

where Di cx £qi and D p cx £opi an d L is the length of the Fabry-Perot resonator of Figure 2. It was 
shown in [9] that parametric amplification at resonance occurs if 

2C7 #% > 1 (8) 
miv a;;: 

where C/ p is the energy stored in the "pump" mode, and the Qi and Q s are the quality factors of the 
"idler" and "signal" modes, respectively. We can express this in terms of the outside laser power of 
the "pump" beam (^outside) by noting the circulating power inside (Vi ns ide) the Fabry-Perot cavity 
for the "pump" mode is to a good approximation 

^inside ^outside (9) 

7T 

where J- is the finesse of the Fabry-Perot resonator, assuming a high reflection coefficient for its 
mirrors. Solving for the outside "pump" intensity gives the threshold condition for parametric 
amplification 

7tc mLuj 2 s 

threshold =27g^; ( 10 ) 

where L is the length of the Fabry-Perot, and c is the speed of light. For estimates of the threshold 
we select a reflection coefficient |l(Jj r = 0.9999 for the two end mirrors, and r2 = 0.5 for the middle 
mirror (M2 of Figure 3). These values yield a spacing between the two peaks in Figure 5 of f s « 10 
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GHz, which also corresponds to a finesse of J- ~ 15700, a Qi ~ 4.5 x 10 8 for a 700 nm wavelength; 
the quality factor of the "signal" is assumed [TTHT3] to be Q s ~ 10 10 , L = 1 cm, and m is assumed 
to be a free mass [H] on the order of m = 2 mg. With these parameters we find a laser power 
threshold of 

^threshold = 530 mW (11) 

It is strongly noted that these estimates are based on the ideal model described above. However, 
it suffices as an order of magnitude estimate and shows that parametric amplification in such a 
scheme may be possible. Furthermore, parametric amplification of acoustic modes has already been 
experimentally observed within Fabry-Perot cavities [15} I16j. 



VI. THE CASE OF AN ALL-MICROWAVE, ALL-SUPERCONDUCTING RESONATOR 

SYSTEM 



Now suppose that one were to replace the optical "Double" Fabry-Perot system in Figure 4 by 
a microwave "Double" all-superconducting Fabry-Perot system, so that the entire system becomes 
an all-microwave, all-superconducting resonator system. Let us also replace the incoming laser 
"pump" beam by an incoming microwave "pump" beam coming in from the right, and let us then 
recalculate the pump threshold power needed for parametric oscillation. For the case of microwave 
resonators having lengths comparable to the microwave wavelength, the cavity finesse J- becomes 
comparable to the quality factor Q p for the "pump" resonator, so that now the estimate for the 
threshold power becomes 

-p I microwave ^ TflL/Ld s , , 

/ threshold | pump — C {LZ) 
^ip Vs 



Since the system is now an all-superconducting system, we shall assume that [TTrfl3] 

Q P ~ Q % ~ Q s « 10 10 (13) 

Again let us assume a mass m of the pellicle mirror to be 2 mg, and a length L of the resonator to 
be 1 cm, which corresponds to a "signal" frequency oj s of 2ir x 15 GHz, and a "pump" frequency 
Up of 2tt x 30 GHz, in the case of a degenerate parametric oscillator for which uj s = uji = u p /2. 



One concludes from (12) that 



nhresholdip ™™ ^ 0.05 /iW (14) 
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Therefore only microwatt-scale "pump" threshold powers are needed for parametric oscillation for 
an all- microwave, all-superconducting resonator system. This is clearly a feasible experiment to 
perform. 

VII. APPENDIX A: GAUSSIAN MODE OF A SC HEMICONFOCAL RESONATOR 
WITH LONGITUDINAL ELECTRIC FIELDS 

Here we examine in detail the properties of the fundamental Gaussian-beam mode of the elec- 
tromagnetic hemiconfocal resonator sketched in Figure 1, following the methodology of the analysis 
of Ince-Gaussian beams introduced in reference |17j . 

In order to satisfy the conducting boundary conditions at the SC surface of the flat mirror 
located at z = in Figure 1, we seek solutions of the vectorized paraxial wave equation [18] which 
result in a longitudinal electric field vector Ei ong satisfying the conducting boundary conditions at 
the surface of this mirror, which in turn leads to a longitudinal component of the electrical force 
on the charges located at the surface of the mirror. 

The analysis begins with the Maxwell's equations for electromagnetic fields in vacuo [19J 



VxE = -- (15) 
<9E 

V x B = +^o^ (16) 

V-B = (17) 

V-E = (18) 

Assuming that all fields have the same complex exponential time dependence exp (— iut) these 
equations become 

V x E = +iuB (19) 

V x B = -io^ £ E (20) 



V-B = 



(21) 



12 



V • E = (22) 
Taking the curl of the first Maxwell equation (19), and using the second Maxwell equation (20), 



one gets 

V x (V x E) = +ioo (V x B) = u 2 fi e B (23) 

Using the vector identity 

V x (V x E) = V (V • E) - V 2 E (24) 

and using the fact that in the vacuum, V • E = (i.e., that there are no charges present within the 
volume of the resonator), one arrives at the Helmholtz equation for the electric field 

V 2 E + w 2 /x e E = V 2 E + k 2 B = (25) 

where 

k = ojyfjji^EQ = u /c (26) 

is the vacuum wavenumber of the EM wave inside the Fabry-Perot resonator shown in Figure 1, 
with c = l/^//x £o being the vacuum speed of light. 

Likewise, one arrives at the Helmholtz equation for the magnetic field 

V 2 B + ^Vo^oB = V 2 B + £; 2 B = (27) 

Note that the Fabry-Perot resonator configuration of Figure 1 has a dominant axis of propaga- 
tion, namely, the z axis. This suggests that a paraxial wave approximation for the z component of 
the electric field might be useful here. 

An important reason for singling out the z component of the electric field for the Helmholtz 
equation is that the boundary conditions at the flat SC mirror require that the longitudinal (or 
normal) component of the electric field does not vanish at the surface of the mirror, but that 
the transverse (or tangential) components of the electric field must vanish at this conducting (or 
superconducting) boundary. This boundary condition singles out a transverse magnetic mode as 
the fundamental mode of the resonator. 

Let us choose the surface of the flat mirror to coincide with the z = plane. Then from 
the above conducting boundary conditions, we expect that the longitudinal component E z is a 
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maximum at z = 0, and also at z = —L at the conducting surface of the curved mirror, but that 
the transverse components of the electric field vanish at these surfaces. We therefore seek solutions 



of the z component of the the electric field from Helmholtz equation (25), i.e., 



V 2 E Z + k 2 E z = (28) 

that satisfy these boundary conditions. From these solutions for E z , we can derive the transverse 
components of the electric and magnetic fields of the modes of the resonator, in a procedure similar 
to finding the fields of the transverse magnetic (TM) modes of a microwave resonator [19j. In 
other words, the conducting boundary conditions for the configuration of mirrors shown in Figure 
1 demand that the z component of the magnetic field must be zero everywhere on the surfaces of 
the mirrors, i.e., 

B z = (29) 
a condition that is also required by the Meissner effect in SC's. 



It is then natural to try as a solution of the Helmholtz equation ( 28 ) the following Ansatz for 
the z component of the electric field: 

E z (x, y, z, t) = Eo^ (x, y, z) cos kz exp (— iojt) + c.c. (30) 

where Eq is a constant, ip (x, y, z) is a dimensionless factor with a complex amplitude that varies 
slowly with z compared to the fast z dependence of the factor cosfcz, and "c.c," as usual, is the 
complex conjugate of the previous term. We have chosen here the z dependence to be that of a 
cosine rather than a sine function in order to satisfy the boundary condition that E z is maximum 
(i.e., an anti-node) at z = 0. This leads to a standing- wave solution for the mode of the SC 
resonator. 

Therefore neglecting the second derivative of if) (x, y, z) with respect to z, which is assumed to 
be small compared to the first derivative with respect to z, we arrive at the paraxial wave equation 
for the slowly varying amplitude 

2ik^ + V,V = (31) 
oz 

where the transverse Laplacian is 
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Note that the paraxial wave equation (31) has the same mathematical form as the time- 
dependent Schrodinger equation 

hdib H 2 o , . 

-d = * < 33) 

except that the axial distance z in the paraxial resonator problem has now been replaced by the 
time t in the quantum mechanics problem. Since we know that spreading Gaussian wavepacket 



solutions are solutions to the time-dependent Schrodinger equation ( 33 ) , we expect that analogous 



Gaussian solutions will also be solutions to the paraxial wave equation (31). 

For the Fabry-Perot resonator configuration of Figure 1, the procedure is to first solve the 



paraxial wave equation (31) for E z after setting B z = everywhere. Then one substitutes these 
solutions into the right-hand sides of Maxwell's equations as source terms to obtain solutions for 
the transverse fields E x and E y . 



The lowest-order Gaussian solution to the paraxial wave equation (31) can be obtained starting 
from the following "Gaussian" Ansatz: 

if) (x, y, z) = * G (x, y, z) = exp i (p(z) + (z) r 2 ^j (34) 

where P (z) and Q (z) are complex functions of the axial distance z, and where 

r 2 = x 2 + y 2 (35) 

is the square of the radial distance r of a field point displaced away from the axis of the resonator. 
Note that r is a real variable. 

It follows from this Gaussian Ansatz that 

^ G = iQ^ G - Q 2 xH G (36) 



I^VEV, - iQVo -(f,,H> ( , (37) 



so that the paraxial wave equation, which is a PDE, reduces to the ODE 

- 2k— - 2k ■ 1 °^r 2 + 2iQ - Q 2 r 2 = (38) 
dz 2 dz 

Collecting the coefficients of r 2 in this equation and setting their sum equal to zero, one then arrives 
at the two first-order ODE's 

k^ + Q 2 = (39) 
dz 

dP 

-k—+iQ = (40) 
dz 
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The solution of (39) can be obtained by integration as follows: 



Transforming variables using 



dQ 1 . , 



(41) 



Q 



(42) 



one finds that (41) becomes 



dq + — I dz = 



</il 



k 



(43) 



and therefore that the solution to the ODE for Q (39) can be rewritten as follows: 



q 00 - qo = 



(44) 



It should be kept in mind that since Q (z) is a complex function of z, so likewise q (z) will also be 
a complex function of z. 

In analogy with the Gaussian wavepacket-spreading problem in quantum mechanics, it is natural 



to impose on the Gaussian Ansatz (34) as an initial condition (i.e., boundary condition) at z = 0, 
that it reduces to the real Gaussian function 



vlV; (•»•.//.()) - exrw | P(0) + ^Q(0)r 2 

= exp (-v 2 /wq) 



(45) 



where the real number wq is the initial Gaussian wavepacket size evaluated at the z = plane, 
which corresponds to the "beam waist" size at the flat mirror in Figure 1. Then 

1 i o 



0(0) 



(46) 



is determined to be an imaginary number. Also, it then follows from the initial condition (45) that 
the initial value of P (z) must be 



P(0) = 



(47) 



Therefore the solution for q (z) MM becomes 



q (z) = qo + 



k 



-2< + k 



z z — ikw\j2 z — izR 



k 



k 



(48) 
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where the "Rayleigh range" zr is defined as follows: 

z R = kw 2 j2 (49) 

The meaning of the Rayleigh range zr is that it is the distance scale on which the width of the 
Gaussian beam will have significantly increased due to diffraction (i.e., wavepacket spreading due 
to the uncertainty principle). The paraxial approximation holds when diffraction angle ^diffraction 
due to the spreading of the Gaussian beam is sufficiently small, i.e., when 

^diffraction ~ — < 1 (50) 
Therefore the paraxial approximation is satisfied when the Rayleigh range satifies the condition 

z R = kw 2 j2 = ttwI/X ~ vrA/^. ffract . on > A/2 (51) 



Transforming back from the function q (z) to the original function Q(z), one finds from (48) 
that 



' (z) \z — izr) \z 2 + z 2 - 



R 



1 , . ZR 



k TT +i ■> 2 ) ( 52 ) 

\z + z 2 R /z z 2 + z 2 R J 



It therefore follows that in the Gaussian Ansatz (34), the exponential factor exp {%Q (z) r 2 ) can 
be factorized as follows: 

( i 2 \ ( % kr 2 

ex P \7;Q( z ) r = ex P 



2~° w ' J ~" r \2z + z 2 R /z 



x ex P ( ~\tT— k -TTY\ I ( 53 ) 



2(1 + ^/4) Z R, 
i kr 2 \ ( r 



exp - exp 



2R(z)J K V w 2 (z) 



,.2 



where 



R(z) = z + z 2 R /z (54) 
is the radius of curvature of the phasefront of (x, y, z) at z , and where 



"•V-!' 1 • ) (•>"») 

Z R 
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is the square of the spreading Gaussian beam width at z. Thus one concludes that the area of the 
Gaussian beam, as measured by w 2 (z) evaluated at z = zr, will double due to diffraction from its 
initial value w 2 . evaluated at the z = plane. 



Next, we shall find the solution of the ODE (40) for P(z), starting from the known solution 



(52) for Q (z), by direct integration, as follows: 



i l' z l' z I 

P( z ) = - Q(z)dz = i dz 

k Jo Jo z- iz R 



i In 



z - iZR 

-izR 



(56) 



Therefore it follows that the first exponential factor in the Gaussian Ansatz (34) can be expressed 
as 



where 



exp (iP (z)) 
exp ^— In 
1 



z 2 + z% 



ZR 



exp (itp G (z)) 



= exp I —i arctan — 
l + z 2 /z 2 R ' V Z R 



(57) 



tfQ (z) = — arctan 



zr 



(58) 



is called the "Gouy phase shift" . 



Putting everything together, one then finds that the Gaussian Ansatz (34) becomes 

>,2 



^c(x,y,z) 



exp 



w (z) *~ \ w 2 (z) 
kr 2 



x exp < i 



- i arctan — 
2R(z) z R 



(59) 



which is a solution of the paraxial wave equation 



2ik 



d^ G d 2 ^ G d 2 V G 



+ 



+ 







(60) 



dz dx 2 dy 2 

and that therefore the full standing-wave solution for the longitudinal electric field in the transverse- 
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magnetic TMoin mode of the SC resonator becomes 

E z (x,y,z,t) = Eq^g cos kz exp (— iut) + c.c. 

TP W ° ( _J^_\ 

w {z) \ w z [z) J 

( kr 2 z 
x exp < i — , . — i arctan — 
y \ 2R{z) z R 

x cos kz exp (— iut) + c.c. (61) 
The first boundary condition at the flat mirror at z = will be satisfied, since 

cos kz = 1 at z = (62) 

and therefore at z = 

f r 2 \ 

E z (x, y, 0, t) = Eq exp k exp (— iujt) + c.c. (63) 

i.e., that the longitudinal electric field has a standing- wave anti-node at z = 0. If one further 
assumes that the curvature of the curved mirror in Figure 1 matches the curvature of the curved 



phase front R(z) given by (54) evaluated at z = —L, then the second boundary condition at the 



curved mirror at z = —L will be satisfied when 

cos kL = ±1 for kL = nir where n = 1, 2, 3, ... (64) 

i.e., that the longitudinal electric field has a standing- wave anti-node at z = —L (we neglect here a 
small correction factor arising from the Gouy phase shift). Therefore the standing- wave anti-node 



condition (64) determines the TMoin eigenmode frequencies of the SC resonator depicted in Figure 



1 through the relationship 



Q 

uj n = ck n = 2-7T • n ^ where n = 1, 2, 3, ... (65) 



Hence the fundamental eigenmode, i.e., the one with the lowest resonance frequency, corresponds 
to the mode number n = 1, i.e., the TMon mode. 

In order to find the transverse components E x (x,y,z,t) and E y (x,y,z,t), let us return to 
Maxwell's equations. Since no free charges are present in the vacuum between the two SC mirrors 
of Figure 1, it follows that V • E = 0, and therefore that 

dE x dE y _ dE z 
dx dy dz 
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From the paraxial approximation condition (51 ), one concludes that the leading term on the right 
hand side will be given by 

dE 



dz 



Eq^/g(x, y, z)k sin kz exp (— iut) + c.c. 



This suggests that we try as trial solutions 



Eq 



E x = x ■ —^ck sin kz exp (—iut) + c.c. 



(67) 



(68) 



E y = y ■ —^ck sin kz exp (—iut) + c.c 



because then 



d(x^ G ) _ d^ G _ 

ox ox 2 



(69) 



(70) 



where the last term can be neglected in the paraxial approximation, so that the leading terms are 

d(x^ G ) 



dx 



(71) 



d(y$ 



G) 



dy 



(72) 



Therefore (66) becomes 



2 dx 2 dy 



(73) 



which holds true in the paraxial approximation. 



As a further check of the validity of the above trial solutions ( 68 ) and ( 69 ) , let us verify that 
Faraday's law in component form 



is satisfied. One finds that 



3E X dE y = dB z = Q 
dy dx dt 



dE x d (x* G ) 

oc ^ ' = x-2y iQ^G 



(74) 



dy dy 



(75) 



(76) 
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cancel each other, so that we see that (74) is indeed satisfied. We therefore conclude that, in the 



paraxial approximation, (68) and (69) are indeed the unique solutions for the transverse components 
E x and E y that correspond to the paraxial TMoi n , Gaussian-beam solution E z given by (61). 



The two transverse magnetic field components B x and B y can then be gotten from the two 
Maxwell equations 



dB x dB,, 



dx 



dy 



dB z 
dz 







dB x dB v dE z 



dy 



dx 



By inspection, the solutions for B x and B y of (77) and (|78j) are 

En 

B, 



U ■ ^-i&Qkcoskz (-ie lujt )+c.c. 



B„ 



2c 

-x ■ — -^ak cos kz (—ie 
2c K 



-iuit 



) + 



C.C. 



(77) 
(78) 

(79) 
(80) 





(b) 



(c) 



FIG. 6: (a) Sketch of a snapshot of the pattern of supercurrents j (black arrows) flowing on the surface 
of the SC flat mirror at z = of Figure 1, and their associated magnetic field vectors B (blue arrows) 
immediately outside of the SC surface. This pattern of currents and fields is associated with the circularly 



symmetric TMoi n mode solution (61). (b) Transverse intensity profile of this mode (i.e., the sum of the 



absolute squares of E x and E y given by (68) and (69), respectively), (c) Longitudinal intensity profile (i.e., 



the absolute square of E z given by (61 )). 



The above solutions can be depicted as the circularly symmetric, transverse magnetic TMoi„ 
mode pattern sketched in Figure 6(a), which is a snapshot taken looking down the z axis of the 
hemiconfocal resonator at the z = plane of the flat SC mirror of Figure 1, at the moment of 
maximum magnetic field. In Figure 6(b), the transverse intensity pattern of this mode is depicted. 
The circular symmetry of this pattern can be readily understood by taking the absolute square of 



the transverse electric fields of this mode, which are given by (68) and (69), i.e., 

E r oc x = r cos 9 



(81) 
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E y oc y = r sin B (82) 

so that 

\E X \ 2 + \E y \ 2 (x r 2 (83) 

which is clearly independent of the azimuthal angle 0. Figure 6(c) shows the longitudinal intensity 
profile of l-E^I 2 . 

VIII. APPENDIX B: SIMPLE HARMONIC OSCILLATOR MODEL FOR THE 
MICROWAVE CAVITY-PELLICLE MIRROR SYSTEM 

Here it is shown that the superconducting (SC) microwave cavity with a pellicle end-mirror (see 
Figures 1 and 4) can be modeled as a simple harmonic oscillator whose loaded quality factor Qioaded 
is approximately given by the quality factor of the SC cavity Q s . 

Let the pellicle end-mirror consist of a thin SC film deposited on a thin, light, flexible diaphragm, 
which is sufficiently thin so that it can easily be driven into mechanical motion. Furthermore, 
suppose that the SC film is electrostatically charged with a net DC charge q. We assume that 
the charge q which resides on the surface of the film is so tightly bound (via the Coulomb force) 
to the metallic film that when the charge q moves, the film will co-move with it [20]. Then the 
longitudinal electric field E 2 at the surface of the SC film will lead to the instantaneous force 

F z (t) = (7E 2 (t) + F rad (t) (84) 

where E 2 is the longitudinal electric field at the surface of the SC film, and where the force on the 
film due to radiation pressure is given by 

F rad (t) = ^e B z (t) 2 A oc E 2 (t) 2 (85) 

where eo is the permittivity of free space, and A is the area of the film over which E 2 (i) is nonva- 
nishing, i.e., the Gaussian-beam- waist area of Figure 6(c). Since the radiation force F ra( j (t) scales 
quadratically with the electric field at the surface of the film, while the Coulomb force qE z (t) scales 
linearly, there exists a maximum electric field strength E mSuX such that if \E Z (t)\ < E mSuX , then the 



Coulomb force qE z (t) dominates over th.6 radiation force Frad 

(t) . Comparing ( |84[ ) and ( [85] ), one 

finds that 

F 

e A 
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Expressing this condition in terms of a maximum externally applied "seed" power (V™*) and the 
maximum circulating power in the cavity (T^v*) we 

« x = ^ « 7 nW (86) 

PSv X = ^ « 43 W (87) 

where we have assumed perfect coupling between the externally applied "seed" power and the 
cavity, A = ttwq is the cross-sectional area given by the beam waist (with wq = 1 cm) of the TMqii 
mode in cavity, q « 20 pC for 100 Volts DC, u = 2vr x 12 GHz, and Q s = 10 10 [12]. For the rest 
of this analysis, we assume that we are in the regime where the circulating power in the cavity is 
sufficiently less than 43 W, or, more generally, that it is less than V^ x , so that the radiation force 



is negligible, and (84) becomes 

F, (t)^qE z (t) 

At high (i.e., microwave) frequencies we take the approximation that the pellicle end- mirror 
behaves like a free mass, so that the equation of motion for the pellicle mirror is given as 

m^=qE(t), (89) 

where we drop the subscript z from ~E z (t) for convenience, and we switch from z to the variable x 
to denote the displacement of the oscillating mass m from equilibrium. The time-dependent part 
of the longitudinal electric field at the surface of the SC film can be described as a harmonically 
time-varying field given by 

E(t) = £{t)e~ wt + c.c. 

where £{t) is a slow varying amplitude. It follows that the displacement of the charged mirror is 
given by 

x(t) = %£(t). (90) 

Observe that the displacement x(t) is linear with the longitudinal electric field E(t) evaluated at 
the surface of the flat mirror in the SC resonator. Now suppose that the SC resonator is in steady 
state and filled with some constant input power from some "seed" microwaves so that the pellicle 



end- mirror displacement is given by (90). If we then shut off the injected "seed" power, the SC 
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resonator's electric field will decay exponentially with time. Hence, the displacement of the pellicle 
end- mirror will also decay exponentially with time. It now suffices to show that the electric field 
E(t) in the resonator can be described by a simple harmonic oscillator. The equation of motion for 
the undriven simple harmonic oscillator is 

d x dx 9 /^-\ 
W + 2 ^t +UJX = ° (91) 



where 7 is the decay parameter of the oscillator. Using (90) we arrive at an equivalent simple 
harmonic motion equation for the field in the cavity 

d 2 E dE 9 

_ + 2 7 -W £ = (92) 

where E is the electric field evaluated at the surface of the moving mirror, where we interpret 
27 as the FWHM of the SC cavity resonance, and where u is the resonance frequency of the SC 



resonator. (Note that (92) also follows from the Helmholtz analysis for a lossy resonator.) In the 



slowly varying amplitude approximation, (92) reduces to a first order linear differential equation 
for the slowly varying amplitude 

f +7 £=0 (93) 

whose solution is 

£{t) = Eoe^ 1 = E e^ t/2Qs (94) 

where Eq is the initial electric field amplitude at the surface of the mirror, and Q s = is the 

SC resonator's intrinsic quality factor. Therefore, the field in the SC microwave resonator decays 
like a simple harmonic oscillator with a time constant that is proportional to the quality factor of 
the cavity. Furthermore, since the displacement of the pellicle end-mirror is linear with the field 



inside the resonator, it must also decay like a simple harmonic oscillator. Note that with (94) we 
can write the field inside the resonator as 

E{t) = £ Q e- w */2Qs e -i^ + cc (95) 
which is the well-known exponentially decaying solution with the ringdown time 

r r = 2Q s /uj 

of the resonator. 



24 



Finally we take into consideration the effect that driving the pellicle end-mirror with "seed" 
radiation has on the quality factor Q s of the SC resonator. The loaded quality factor Qi aded> 
where the loading refers to the power loss due to the simple harmonic motion of the charged 
mirror, is given by 

Qloaded = - (96) 
' loss 

where Uq is the energy stored in the cavity, and 'Pioss is the total power loss in the cavity given by 

^loss — T^c ^mirror (97) 

where V c is the intrinsic power loss of the SC resonator and is related to the resonator's intrinsic 
quality factor by Q s = uiUq/V c , and where "Pmirror is the average power loss due to the motion of 



the charged pellicle end- mirror. From (96) one finds that 



n _ QsQshm (Q „, 

(^loaded — T7S V y8 J 

+ VSHM 

where Qskm = ^tW-Pmirror is the contribution to the quality factor arising from simple harmonic 
motion (SHM). Although it is possible that some or all of the power loss that goes into the simple 
harmonic motion of the charged mirror is converted into electromagnetic radiation power which 
goes back into the SC resonator, it is instructive to account for it. The average power loss due to 
the moving pellicle end-mirror is 

^mirror = (F • v) = (99) 

The electric field is calculated by assuming some externally applied "seed" power V ey ± is injected 
into the SC resonator. In steady state, the energy in the cavity L^o is 

tt _ jffijgxt Qs Mnm 
P °-(l+fV (100) 

where f3 is a coupling parameter of the input/output hole in Figure 1 and is assumed to be unity, 
(3 = 1. It follows that the amplitude of the electric field inside the cavity E$ is given by 

2 = (1Q1) 

where V is the effective Gaussian-beam volume of the SC resonator. Hence, the average power loss 
from the pellicle end-mirror is 

v _ kfVextQa (W2) 

' mirror — 9, , V ±u-i / 

mu z Veo 
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Assuming an external applied "seed" power of V e *t = 500 pW, q = 20 pC, Q s = 10 10 , m = 2 
mg, oj = 2tt x 12 GHz, and an effective Gaussian-beam volume of V = 1.25 cm 3 , one finds 

^mirror » 6 X 10~ 2 ° W (103) 

This power yields 

Qshm ~ 4 x 10 21 . (104) 



Since Q s <C Qshm , it follows from ( 98 ) that to an extremely good approximation 

Q loaded ~ Qs 

IX. APPENDIX C: SIMPLE SOLUTION TO A FABRY-PEROT WITH A 
HARMONICALLY MOVING END-MIRROR IN QUASI STEADY STATE 

We consider a Fabry-Perot under quasi steady state conditions and ignore transients during the 
build up of the modes as those discussed in [21]. In a quasi steady state, the harmonically moving 
end- mirror generates two Doppler side bands, "Stokes" and "anti-Stokes" |22| . The "anti-Stokes" 
sideband is suppressed via the "Double" Fabry-Perot scheme as illustrated in Figures 3 and 5. 

The radiation force is given by 

F(t) = ^e \E\ 2 A (105) 

where eo is the permitivity of free space, A is the cross sectional area of the Gaussian beam, and E 
is the total electric field in the Fabry-Perot. Modeling the moving end-mirror as a simple harmonic 
oscillator as depicted in Figure 2, we have 

•> Fit) , , 

x + 2j s x + u} 2 s x = — 106 
m 

where x is the displacement of the simple harmonic oscillator from equilibrium, uj s is the natural 
oscillator frequency, m is its mass, and 2~f s is the FWHM. Using the slowly varying amplitude 
approximation in quasi steady state 

x = X{t)e~ iWst + c.c. (107) 



the left side of (106) becomes 



2/>„ ( d ^- + j s X) e~ iWst + C-.C-. ( .1.08 ) 
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where we have assumed that 2 (j s — ico s ) ~ —2w s - The right hand side can be expanded in terms 
of the fields 

E i =£ 0i {t)e- i ^ t + c.c, E p = £ 0p e- i "P t + c.c. (109) 

where the former is the "Stokes" or "idler" term and the latter is the "pump" mode, £oi(t) being a 
slowly varying amplitude in quasi steady state, but £q p being a constant, in the "undepleted pump" 
approximation. Taking the beat terms and neglecting the nonresonant terms \Ei\ 2 and \E p \ 2 we 
have 

F(t) =e A[r 0i £ 0p e- i ^-^ t 

+ SoiSope- 1 ^^ 1 + c.c] (110) 

Equating both sides and multiplying by e luJst we find 

-2iu s ( d ^ + ls x\ +c.c.(oc e 2 ^) 



m 

+c.c.(oce i2uist ,e i2 ^% (111) 



where uj p — uoi ~ oj s ■ Here we employ the fact that this is a linear system and observe that only the 
force at resonance will be the main driving force of the mechanical oscillator; hence, we can neglect 
off-resonance terms (in the rotating- wave approximation), and write the driven oscillator equation 

as 



-^^^ = 2^°^ " " " (112) 
Following |9], let us define the fields as 

E i (t)=A i (Di(t)e- iUit + c.c.) (113) 

E p = ApiDpe-^ + c.c.) (114) 

where Ai, A p are normalized so that the total energy stored in each mode is U Pt i = 2uj P: i\D P: i\ 2 , 

and Di {t) is a slowly varying complex amplitude, but D p is a constant, in the "undepleted pump" 
approximation. Calculating the energies in each mode we find 

U p ,i = ^e (E 2 ) = A 2 pi eo\D Pt i\ 2 AL = 2u 2 p ^D p tf (115) 
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where V = AL is the volume of the cavity. Solving for A p; , 



e V 

Inserting the normalization 



A P,i = U1 p,i\l (116) 



£* Qi = AiD* and £ 0p = A P D P (117) 



into equation (112) we arrive at the result [9], eq. (2), 

— +ls X = —5- 
at mu s L 

where Auj = uj p — Ui — uj s . 



+ lgX = l -^i Dp D*e- iA ^ (118) 



To arrive at the equation for the "Stokes" or "idler" field in the cavity we note that the "Stokes" 
mode is generated from the main "pump" mode so that the "Stokes" field is proportional to the 
oscillator amplitude. Recall that we are assuming quasi steady state conditions. The electric field 
reflected from a moving mirror is given by [22] 

E p = £ 0p e- luJ ^e 2ik ^ x + c.c. (119) 

For small x, we find 

E p = £ 0p e~^*(l + 2ik p x) + c.c. (120) 

where the last term leads to the Doppler-generated electric field 

^Doppier = 2ik p x£ 0p e- tuJ " t + c.c. (121) 

Since we assume quasi steady state conditions, the "pump" optical mode is the main source of the 
"Stokes" optical mode. As is well known, the Fabry-Perot cavity field obeys a recursion relation; 
for example see [2H [23] . Hence 

Ei(t + t) = E Dovplcr (t + r) + KP f Ei(t) (122) 

where r = 2L/c is the round trip time, 1Z is the power reflectivity (i.e., the absolute square of the 
reflection coefficient) of the end mirrors, Pf is a propagation factor |21| [23] which accounts for the 
phase accumulated by the beam one round trip earlier, and -Eboppier is the Doppler electric field 
generated from the "pump" mode. Putting in the slowly varying displacements and fields as before 

x = X (t) e~ luJst + c.c. (123) 
Ei = £ i (t) e"*" 1 * + c.c. (124) 
E p = £ 0p e- iuJpt + c.c. (125) 
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and taking the resonant terms, we find 

d ^ +li £ 0l = l -X*£ Qp e-^ t (126) 
where we make use of Aw = uj p — Ui — oo s . Inserting the normalization used in [9] we find 

^ +7iA = L^ X ^ e (12?) 

where we have used the fact that r = 2L/c, and k p c = lj p . With the approximation that 

cj? cj? u;H 

= 5_ ~ _p = Wp (128) 

— w s w p 

(i.e. w p > w s ) we find, 

d ^ + li D i = iX *y* e-^. (129) 
This result is consistent with [Sj, eq. (1). For the detailed solutions to the coupled differential 



equations (118) and (129) and the condition for parametric amplification the reader is referred to 



El. The threshold condition is 

2UpQiQ s > l 



mL 2 uj^ 

where U p is the energy stored in the "pump" mode, L is the length of the Fabry-Perot, m is the 
effective mass, and the Q's are the quality factors. 
Note that the Q's are defined as follows 

Q s = — , "Mechanical" Oscillator (131) 
Q; = — , "Stokes" mode in FP (132) 
Q p = "Pump" mode in FP (133) 

where the 7's correspond to the relaxation rate for each of the modes (i.e HWHM). 
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